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Abstract 

A non-isotropic version of phase equations such as the Burgers equation, 
the K-dV-Burgers equation, the Kuramoto-Sivashinsky equation and the Ben- 
ney equation in the three-dimensional space is systematically derived from 
a general reaction-diffusion system by means of the renormalization group 
method. 

PACS codes: 47.20.Ky 

Keywords: phase equation, renormarization group method 



1 Introduction 



It was a long time ago that the Kurainoto-Sivashinsky (K-S) equation was 
proposed as a higher-order phase equation describing an unstable phase state 
|jl[]. Although there have been many works based on the K-S equation since 
then, there are only few derivations of the K-S equation based on the singu- 
lar perturbation methods. Recently, the isotropic K-S equation has been de- 
rived as a phase equation of a periodically oscillating solution of the complex 
Ginzburg-Landau equation by means of the renormalization group method 
0]. In a reaction-diffusion system, spatial symmetry often breaks so that 
both spatially and temporally oscillating solutions emerge. In such a case, it 
is anticipated that a slowly varing phase of the symmetry-breaking oscillat- 
ing state is asymptotically governed by a non-isotropic phase equation. The 
Gross- Newell's phase equation was derived by the RG method 0. However, 
there are no explicit derivations of such non-isotropic phase equations as the 
non-isotropic K-S equation, the K-dV-Burgers equation and the Benney equa- 
tion 1^ from a general reaction-diffusion equation. In this paper, we derive 
such non-isotropic phase equations to a symmetry-breaking state of a gen- 
eral react ion- diffusion equation by means of the renormalization group (RG) 
method. 



2 Renormalization Group Method 



The perturbative RG method introduced in [Q is shown to be interpreted as 
the procedure to obtain an asymptotic expression of a generator of a renormal- 
ization transformation based on the Lie group ^. This Lie group approach 
provides the following simple recipe for obtaining an asymptotic form of a RG 
equation from ordinary differential equations (ODE). 

(1) Get a secular series solution of a perturbed equation by means of naive 
perturbation calculations. 

(2) Find integral constants, which are renormalized to elliminate all the secular 
terms in the perturbed solution and give a renormalization transformation for 
the integral constants. 

(3) Rewrite the renormalization transformation by excuting an arbitrary shift 
operation on the independent variable:t — > t + r and derive a representation 
of a Lie group underlying the renormalization transformation. 
(4) By differentiating the representation of the Lie group with respect to arbi- 
trary r, we obtain an asymptotic expression of the generator, which yields an 



asymptotic RG equation. 

This procedure is valid for general ODE regardless of translational symmetry. 
The above recipe for ODE is also applicable to autonomous partial differn- 
tial equations (PDE) by choosing suitable polynomial kernels of the linearized 
operator appearing in perturbed equations. First, we should take the lowest- 
order polynomial, of which degree is one, as the leading order secular term. As 
perturbation calculations proceed to the higher order, polynomial kernels of 
higher degrees are included in the higher-order secular terms order by order. 
When we can continue this process consistently, we say that secular terms 
of polynomials are renormalizable or , simply, the consistent renormalization 
condition is satisfied in the sense of the Lie approach. If this is the case, we 
can determine suitable polynomial kernels among infinite number of kernels of 
the linearized operator and the step (1) in the recipe is completed. There are 
no problems in the other steps. Thus, the Lie-group approach is consistently 
applicable to PDE and derivation of some soliton equations and simple phase 
equations (e.g. the isotropic Burgers equation) was presented in 0. 
In the following sections, using this Lie-group approach of the RG method, we 
derive various non-isotropic phase equations from a general reaction diffusion 
system. 



3 Non-isotropic Burgers Equation 

Let us consider a general reaction-diffusion system of equations: 

dtU = F{U) + DV'^U, (3.1) 

where U is an n-dimensional vector and D is an nxn constant matrix. Suppose 
( ^.1| ) has a spatially and temporally oscillating solution Uq = UQ{k, kx—uj{k)t+ 
^Jsatisfying 

- iujUo,e + F{Uo) + k^DUoM) = 0, (3.2) 

where k and are arbitrary constants; uj{k) is a definite function oi k; 6 = 
kx — ut + 0; the suffix 9 denotes the derivative with respect to 9 and Uq is a 
2tt periodic function of 9. For later convenience, we list some useful identities. 
Differentiating (|3.2|) with respect to 9 and k, we have 



L{9)Uo,e = -{u;de + F'{Uo)-+eDd's)Uo,e = 0, (3.3) 

LUo,ee = F" : Ul^, (3.4) 

Lf/o,fc = (u + 2kDde)Uo,e = MUo,9, (3.5) 

LUo,kk = F" : (f/o,fc)2 + 2Mf/o,efc + 2DU^^e0 + ^U^,e, (3.6) 

LUofik = MUo,ee + F" : Uo,eUo,k, (3.7) 



where F' -V = {V ■Vu)F{U)\u=u,, F" : VW = [iW-\/u)iV ■Vu)F{U)]\u=Uo, 



uj = dkOJ and u = dlu. 



Let us seek a secular solution close to U^ik, 6): 

U = Uo{k + k{x, r^, t),9 + 6{x, r^, t)) + U{9, x, r^, t), (3.8) 

where r± = {0,y,z); 6{x,r±,t) and K{x,rj_,t) are small secular deviations 
from the constant phase (p and the wavenumber k respectively and so 

dx6 = 6x = K,. (3.9) 

lf{6,x,rj_,t) represents small perturbed fields which modifies the 0-th order 
field pattern Uq and is not expressed by differentials of Uo{k, 6). The arguments 
(a;, r_L,t) designate "secular variables", that is, all secular perturbed fields are 
polynomials with respect to (x, r_|_,t) and periodic with 9. An expansion of 
f/o in terms of small deviations 5 and n yields 

Uo{k + K,e + 5) = Uo{k, 9) + 5Uo,e + KUo,k 

+ {6y2)Uo,ee + {^V2)Uo,kk + ^^U^fik + ■■■■ 

Note that 5{x,r^,t) (and k) is a secular or poljTiomial function of (a;, r_|_,t), 
which should be eliminated by renormalizing the phase later. 
In this section we suppose that 5 and U are expanded in terms of a small 
perturbation parameter e as 



5 = 


= e(Pi + eP2 + ...), «: = 


= e(Pi,, + eP2,x + ■ 


■■) 


(3.10) 


U = 


= e'(t/2(^) + ef/3 + ■■■), 






(3.11) 



where the suffix x denotes the derivative with respect to x and Pj {j = 
1,2,---) are polynomials of {x,rj_,t), which have increasing degrees with j 
so that polynomial secular terms are renomalizable in the sense of the Lie 
approach of the RG method 0. U2 depends only on 9 since the leading order 
term of U is not contain secular terms of {x, r±,t) ; otherwise they would not 
be eliminated by the RG procedure (see ( |4.58|) and (^4.59|) or ( ^.90|) and ( p.91| ) 
). Then (|3.8|) reads 

U = Uo{k,9) + eUi + e''U2 + ---, (3.12) 

Ui = PiUo,e{k,d) + Pi,.UoAk,0), (3.13) 



U2 = P2Uo,e + P2,.Uo,k + {l/2)PtU, 



0,09 



+PiPi,.Uo,ek + {l/2)PlJJo,kk + U2. (3.14) 

Introducing the Galilean transformation 

x = X — ut, t' = t, (3.15) 



and substituting ( |3.12| ) into (pTT|), we have to the first order perturbed field f/i 
{dt> + L{e) - 4,M - V'^D)Ui = LUi = 0, (3.16) 

where {6,x',t',r±) is considered as a set of independent variables. Here after, 
the prime attached to {x', t') is omitted for simplicity. Substituting (|3.13| ) into 
( P.16| ) and using ( p.3|) and (|3.5| ), we have 

+Pi,tf/o,0 + Pi,,tf/o,fc = 0. (3.17) 

Since Pi is a polynomial of (x, r±, t) and Uq is a periodic function of 6 , ( p.l7| ) 
reads 



l,a;a; 



0, V'Pi=0, V^Pi,, = 0, Pi,t = 0, Pi,.t = 0, (3.18) 



where the suffix t denotes the derivative with respect to t. Noting that the 
leading order secular term consists of a polynomial of degree one, (|3.18|) yields 



Pi = Pi,.x + VxPi ■ rx. (3.19) 

where V_l = (0,(9y,(9z); Pi^x and V_lPi are arbitrary constants. 
The second order equation obeys 

~LU2 = (1/2)P" : Ul{e), (3.20) 

Substituting ( ^141) into (lOOl) , we obtain with the aid of (^-(^ 



+V^P2DUo,ei0) - P2,tUo,eid) " ^2,.,tf/o,fc(^), (3.21) 

which is an equation for periodic U2- Since U2 is a function of 6* only, all the 
coefficients of functions of 6 in the right hand side (RHS) of ( |3.21D do not 
depend on (x, r_L,t), that is, 

P2,xx = Cl, V P2 = C2, P2,t = C3, P2,xt = C4; (3.22) 

where ci,C2 and C3 are non-zero constants while C4 = due to the consistent 
renormalization condition, i.e. both t and xt do not enter in P2 as secular 
terms to be removed consistently by renormalization. This requirement holds 
throught this paper. An explict form of P2 is given as a polynomial of degree 
two with respect to (x, r_|_): 

P2 = P2,xxxy2 + {r^r^ : V^ Vx)P2/2 + x(r^ : Vx)P2,. + Ps.^t, (3.23) 



where all the coefficients of monomials, i.e., -P2,aa-, V_lV_lP2, V_LP2,a: and P2,t 
are arbitrary constants and 

" " k=n / \ 

(fTTTTTT : V^ ■ ■ ■ Vx) = E U ) ?/"~'^'^r'^'- (3-24) 

Then, a periodic solution U2 is possible only if the following compatibility 
condition is satisfied. 

\V±Pi\^{U ■ DUo,ee) - (l/2)u)(Pi,,)2(?7 • Uo,e) + P2MU ■ MUo,k) 

+V2p2(f/ • DUo,e) - P2AU ■ Uo,e) = 0, (3.25) 



where U is an adjoint function of a null eigenfunction of L and < U ■ U >i 
A secular solution up to O(e^) is 



+ e(Pi,, + eP2,,)Uo,k + il/2yiP,^^yUo,kk 
+ e2PiPi..f/o 



U = Uo{k,kx-ujt + (f)) + e{Pi + eP2)Uo,e + {l/2)e'PtUo,, 

,x) Uo^kk 
'Ofik 

= Uoik + 4)^,kx-ujt + ^), (3.26) 

where is a renormalized phase defined by a renormalization transformation 

0(x, r^, t) = + S{x, rx, t) = + e(Pi + ePa + ■ • ■)• (3.27) 

Since the renormalized phase should enjoy translational symmery with respect 
to independent variables, ( |3.27| ) is rewritten as , shifting the origin (x, r±,t) = 
(0, 0, 0) to an arbitrary point (x, r_|_, t), 

0(t + r,x + ^,r_L + ?7 ) =0(a;, r_L,t) 

+ e{Pi(e, 1; X, rx, t) + eP2(e, r], r; x, r^, t)X,3.28) 

where (x, r_|_,t) in polynomials Pi and P2 are replaced by {^,t],t) and their 
coefficients depend on the coordinate of the origin (x, rj_,t) ,e.g. 

Pi{^,'n;x,r±,t) = Pi,x(a;, r_L,t)^ + V±Pi(x, r_L,t) ■ r? 

and so on. This reinterpretation of coefficients of secular terms is the key 
ingredient of the Lie approach of the RG method M. Hence, Eq.( |3.28| ) with 
( CT ) and ( p:^ reads 

4>. = [ag(6Pi + e'P2)]o = ePi,., 
Vx0 = [9r7(ePi + e'P2)]o = eVxPi, 
4>xx = [9|(e2p2)]o = e'P2,xx, 

Vi0 = [a^(62P2)]o = 6VIP2, 

0t = [9.(e'P2)]o = e'P2,t, (3.29) 



where [/(^, 77, r; x, rj_, t)]o = /(O, 0, 0; x, r_L, t). Substituing these relations be- 
tween differentials of the renormalized phase and the reinterpreted coeffi- 
cients of polynomials Pi and P2 into (|3.25| ), we obtain a non-isotropic Burgers 
(n-Burgers) equation: 

4>t = 1^110.. + D^Vl4> + iV|i(0x)' + iV^lVx^r, (3.30) 

and 

D\i = D^ + Dl 

D^ = {U-DUo,e)/{U-UoM), D\^ = {U ■ MUo,k) / {& ■ Uo,e) 

N^ = {U.DUo,ee)/{U-Uo,e) 

iV|| = N^ + {{l/2){U-F":iUo,ky) + {U-MUo,ek)}/{U-Uo,e), 

= -u/2, (3.31) 

where the last equality of ( p.31j ) comes from ( p.6|) . 



4 K-dV-Burgers Equation 



In this section, it is assumed that the diffusion coefficient D\\ along the x 
direction in the n-Burgers equation (p.30|) is small as e: 



D|l oc {U ■ iDUo,e + Mf/o,fe)) ~ 0(e). (4.32) 

Nevertheless the net diffusion in the x direction is supposed to be much greater 
than that in the r± direction, i.e. 

Vjd, ~ 0(e). (4.33) 

Eq.( [4.32|) implies that there is a periodic vector V{9) such that 

LV{e) = DUo^e + MUo,k + 0(e). (4.34) 

Suppose that 6, U, and U are expanded as 

6 = e2{Pi(a;,rx,t) + eP2(x,r^,t) + ■■■}, (4.35) 

U = e'{U2{e)+eU3{e,x,r^,t) + ■■■}, (4.36) 

U = Uoik,e) + e^iUi + eU2 + e^U3 + ---), (4.37) 
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then the leading order perturbed terms (0(e^)) gives the same equations as 
( p.l6| )-( ]3.18| ). As a polynomial solution of ( p.l8| ), we choose Pi as 



Pi = Pi^xX, Pi^x = constant, (4.38) 



instead of (|3.19| ) since |-Pi,x| ^ |V_lPi| due to the assumption ( [4.331) . 
To O(e^), we have 

LU2 = 0, (4.39) 

U2 = P2{x,r^,t)UoM + P2A^,'rL,t)UoM + U2{0), (4.40) 



Substituting (|^ into (|]3|), we get 



LU2 = P2,..(M[/o,fe + DUo,e) + VlP2DUo,9 + V^P2,ccDUo,k 

-P2,tUo,e - P2,.tUo,k. (4.41) 

Due to the consistent renormalization condition described in section 2 (i.e., 
consistent increasing of degrees of polynomial secular terms P,), P2 is a poly- 
nomial of degree two with respect to x and of degree one with respect to r^ 
at most so that we can set V\P2 = V^P2,x = and P2^xt = 0, while 

P2,xx = constant = R2. (4.42) 



If P2^t 7^ 0, the fourth term of the RHS of ( |4.41| ) is only a term which causes a 



secular behaviour of U2{0) with respect to 6 and so P2^t = 0. Then by virtue 
of ( [4.34 ) we have 

1/2(0) = R2V{e), (4.43) 

and 

P2 = R2X^/2 + VxP2 ■ r±, VxP2 = constant. (4.44) 

To 0(e'^) we have 

LU3 = (1/2)P" : (t/i)2, (4.45) 

Us = P3{x,r^,t)Uo,e + P3A^.r^,t)Uo,k + {l/2)iPifUo,e9 
+PiPi,.Uo,ek + il/2)iP,,xfUo,kk 
+Us{9,x,r^,t), (4.46) 

from which we get 

LUs = Ps,x.MUo,k + '^^P3DUo,e-{l/2)iP,,xfOUo,e 

+V^Ps,xDUo,k - P3,tUo,e - P3,xtUo,k. (4.47) 



In view of ( [4.38| ) and ( [4.44| ), the consistent renormalization condition requires 
that P3 contains x^ and the first term of the RHS of ( [4.46| ) is secular with 
respect to x. Therefore we set 

U3{e,x,r^,t) = R,{x)V{d) + U;{e), (4.48) 

where 

-R3 = Rs^xX, Rs^x = constant, (4.49) 

and ( |4.47| ) is rewritten as 

LUsie) = -{Rs-P3,x.){MUo,k + DUo,e) + VlP3DUo,e-{l/2){Pi,x)^uUo,e 
-R3,xMV{e) + V^P3,xDUo,k - P3,tUo,e - P3,.tUo,k- (4.50) 

Since ( [4.50|) is an equation to U^IO), all the coefficients of functions of 6 in the 
RHS of ( [4.50|) do not depend on (x, r_|_, t),that is, 

-R3 ~ P3,xx = Ci, Vj^Ps = C2, 

V'P3,x = C3, P3,t = C4, P3,xt = C5, (4.51) 

where c„ are arbitrary constants. The consistent renormalization condition 
yields 02 = 05 = 0, C3 = Ps^xxx and Ci = so that a secular coefficient R3{x) 
is consistently removed by a renormalization transformation to (see ( |4.59| )). 
Thus (^^511) gives 

P3 = P3,xxxxy3\ + x(rx : Vx)P3 + P3A (4-52) 

The compatibility condition for a periodic solution UsiO) requires that 

P2,xx{{U ■ MUo,k) + ■ DUo,e))/e - {u/2){P,^x?{U ■ Uo,e) 
+P3,xxxi{U ■ MV) + {U ■ DUo,k)) - P3AU ■ Uo,e) = 0. (4.53) 

Here the first term in the LHS of ([4.53|) comes from O(e^) terms. 
Now we arrive at a renormalization transformation to O(e^) 

4>ix, r^, t) = + e^{Piix) + eP2(x, r^) + e^P^ix, r^, t)}. (4.54) 

Following the same procedure as in the case of the n-Burgers equation, we 
obtain from (|3§,(|;3B|),(|^) and (|3g) 

(px = ^ Pl,xi 4>xx = e P2,xx, 4>xxx = e P3,xxxj 4>t = ^ P3,t- (4.55) 

Substituting ( [4.55| ) into ( [4. 531 ), we arrive at the K-dV-Burgers (K-B) equation: 

0i = D\\4)xx + A^ii {4>x)^ + A4)xxx, (4.56) 



where 

A = B + A', 

B = {U-MV)/{U-Uo,e), A' = {U ■ DUo,k) / {U ■ Uo,e) ■ (4.57) 

Notice that secular terms in U are also eliminated by the present renormaliza- 
tion procedure. To 0(e''), ( [4.43| ) and ( |4.48| ) give 



U = e'{R2 + eRsix))}V{9) + e^Us{9), (4.58) 

where a secular term R3 is removed by introducing a renormalized R2 such 
that 

R2{x, rx, t) = R2 + eRsix). (4.59) 

This renormalization transformation is found to be consistent with that to the 
phase (|3^, since ( ^^ , (|3|) and (^311) with ci = imply R2 = (pxx/e^- 



If we introduce the following expansion 

6 = e^{P^ + eP2 + ■■■}, 
U = e'{U2 + eUs + ---}, 
U = Uo{k,6) + e\Ui + eU2 + e^U3 + 



instead of ( |4.35| ),( ^75BD ,( [4.37D , and an auxiliary condition D\\ ~ O(e^) instead 



of ( |4.32| ), the similar procedure as above yields , to 0(e*), the following K-dV- 



Burgers equation with the perpendicular diffusion term 

4>t = D\\4)^^ + DxVi0 + N\\ (4)2 + A0,,,, (4.60) 



5 Non-isotropic Kuramoto-Sivashinsky Equa- 
tion 



The diffusion coefficient Dx_ in the n-Burgers equation ( |3.30D is assumed to 
be small so that 

Dxoc(t/-Dt/o,e)~0(e2). (5.61) 

and net diffusions along both x and r±_ directions are nevertheless the same 
order of magnitude, that is, 

dl/Vl ~ 0{e% (5.62) 



Eq.( [5.61| ) guarantees existence of a periodic vector W{6) such that 

LW{e) = DUo^g + 0(e2). (5.63) 

Suppose that 6, U, and U are expanded as 

5 = e'{Piix,r^,t) + eP2{x,r^,t) + ■■■}, (5.64) 

U = e''{U2{9)+eU3{9,x,r^,t) + ■■■}, (5.65) 

U = Uo{k,e) + e^iUi + eU2 + e^U3 + ---), (5.66) 



then perturbed equations up to O(e^) give the same equations as (|4.39| ) and 
(iDfort/, (j = 1,2,3). 

LUj = 0, (5.67) 

U, = P,ix,r^,t)Uo,eid) + PU^,r^,t)Uo4e) + U„ (5.68) 



Noting Ui = and following (|3:T6|), (|3T7|) , and (^J^) , we have (^J^) with 
Pi^x = due to the assumption ( |5.62| ) ,that is, 

Ui=Piir^)Uo,e, Pi = VxPi-rx. (5.69) 

For j = 2, substituting (|5.68|) into (p.67|), we have 



LU^ie) = VlP2DUo,e + P2MMUo,k + DUo,0) + V^P2,.DUo,k 

-P2,tUo,e - P2,xtUo,k, (5.70) 

Since U2 is a function of 6 only, all the coefficients of functions of 6 in the 
RHS of (|5.7(]| ) must be constant. The consistent renormalization condition 
and ( ^.691) irnplies that a degee of polynomial P2 should be two with respect 
to r± and one with respect to x and t at most, i.e. 

P2,xx = V P2,x = P2,xt = 0. 



The first term of the RHS of (|5.70|) does not cause a secular behaviour of U2 



owing to (|5.63| ) and so P2,t = is necessary for a perodic solution f/2. Thus 
we have 

U^{e) = S2W{e), S2 = VIP2 = constant, (5.71) 

P2 = P2,xX + (r^rx : Vx Vx)P2/2. (5.72) 



where all the coefficients of monomials in (p.72| ) are arbitrary constants. 

For j = 3, we have the same equation as ( [5.70| ),where the suffix 2 is replaced 

by 3. 

LtJsie) = VlP3DUo,e + P3MMUo,k + DUo,e) + V^P3,xDUo,k 

-P3,tUo,e - P3,xtUo,k, (5.73) 
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The consistent renormalization condition implies that a degree of polynomial 
P3 should be three with respect to r_L and so V^^Ps in the first term of the 
RHS of ( p. 731 ) is secular with respect to r±, while 

0. (5.74) 







P3,xcc = 0, 


V^Ps,. = 0, Ps,.t 


Therefore, 


U3 


takes the form 








U3 = 


Ss{r^)w{e) + Us{e), 



(5.75) 
where S{rj_) is a polynomial of degree one with respect to r±. Then, ( p. 731 ) 

mie) = (ViPa - S3)DUo,e - P3,tUo,e. (5.76) 

Here P^t = is again necessary for a periodic solution U^^O) and V^Pa — Sa = 
so that a secular term 5*3 is automatically elliminated as soon as the phase is 
r^normalized as shown in the last paragraph of this section. Now, we set 
UsiO) = without loss of generality and have 

^3 = VIP3 = Vx^3 ■ r-x, (5.77) 

P3 = a;(rx : Vx)P3,. + (r^rxrx : VxVxVx)P3/3!, (5.78) 

where V_lS'3, V_LP3,a; and V_l V_lV_lP3 are arbitrary constants. 
A nonlinear term enters in the perturbed equation to O(e^): 

Lf/4 = (1/2)P" : f/2, (5.79) 

f/4 = Piix, rx, t)Uo,e + P^A^, rx, t)Uo,k + il/2)P^{r^)Uo,ee 

+ f/4(^,x,rx,t), (5.80) 

from which, we get 

Lf/4 = V^P^DUo,e + P4,..MUo,k + V^P4,.DUo,k 

+ il/2)V''P^DUo,ee - Pi,tUo,e - P4,.tUo,k. (5.81) 

The similar discussion leading to ( p.75| ) yields 

Ui{0, X, rx, t) = Si{x, r^)W{9) + U^iO), (5.82) 

where 

^4 = 54,.x + {r^r^ : V±V±)Si/2, (5.83) 

and S^^x and V_lV_lS'4 are arbitrary constants. Then (|5.81|) is rewritten as 

LU^iO) = -{S^-V^P^)DUo,e + PA,..MUo,k + S^,xMW 

+V^P,,xDUo,k + {l/2)V^P^DUo,ee + V^S.DW 
—Pi,tUofi — Pi^xtUo^k, 
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which is an equation for [/^{O) and all the coefficients of functions of 9 should 
be constant. 

5*4 — V P4 = Ci, P4:,xx = C2, V P^^x = C2, PA,t = Ca, (5.84) 

and C5 = P^^xt = holds again. The compatibility condition for U4{9) requires 

{V^P2/e^){U ■ DUo,e) + P^MU ■ MUo,k) + VlP,,x{U ■ MW{e)) 
+ViP4,.(f> • DUo,k) + \^±Pi\^{U ■ DUo,ee) + "^IpS ■ DWi9)) 
-P4AU-Uo,e) = 0. (5.85) 

A renormaliaztion transformation to O(e^) is 

0(x, rx, t)=(f) + e=^{Pi + eP2 + e'Ps + e^Pi}, (5.86) 

which gives, noting (|5:69|),(^72D,(^78D, and (p|), 



Vx0 = e'V^Pu 4. = e6P4.,, Vi0 = eViP2 



eV'P 



4,X; 



e^VlP 



4, n 



e'P.f 



(5.87) 



Substituting (|5.87| ) into (|5.85| ), we arrive at a non-isotropic Kuramoto-Sivashinsky 
(n-K-S) equation: 



D'.4>xx + D^Vl4> + iV±| Vx0r + EVi4> + GVl4>x, 



(5.88) 



where 



E = {U-DW{e))/{U-Uo,e), 

G = A' + H 

A' = {U-DUo,k)/{U-Uo,e), H = {U ■ MW (6)) / {U ■ Uo,e) ■ (5.89) 



Secular terms in U are also eliminated by the present renormalization proce- 
dure. To 0{e^), ( pTTD , ( |5T5| ), and ( p:82|) give 



[/ = e4{S2 + eS.ir^) + e^S.ix, r^)}W{9) + 6^,(9), 



(5.90) 



where secular terms 5*3 and 5*4 are removed by introducing a renormalized S2 

m 1 r* h 1 h /^ T 

S2{x, r^, t) = S2 + eSsir^) + e^S^ix, r^). (5.91) 

This renormalization transformation is found to be consistent with that to 
the phase_(/) j ^M ), since (|5Tl| ), ( |5T7D , (|]8^, and (|]8§ with ci = imply 

^2 = Vi0/6^ 
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6 Benney Equation in Three Dimension 



In addition to the assumption ( p.61D , the wave number k is also assumed to 
be as small as 0(e). Furthermore it may be reasonable to assume that uj{k) 
and Uo{k, 6) are functions of /c^, which is satisfied in the case of the complex 
Ginzburg-Landau equation analyzed in the next section. Then the following 
estimates hold 

u^O{e), f/o,fc~0(e), (6.92) 

and 

A'oc (?7-Dt/o,fc) ~0(e), B(x{U ■MV{e)) ^0{t), (6.93) 
D(|Oc(t/-Mt/o,fc)~0(e2), (6.94) 

or 

A = A' + 5- 0(e), D|| =D^ + D|| ~0(e2), (6.95) 

namely both coefficients of diffusion Z^x and D\\ are as small as O(e^) while 

the coefficient of dispersion A is 0(e). Here, the linearized operator L and L 
are rearranged as 

L = L' - k'^Ddl (6.96) 

L' = -Lude-F'{Uo); (6.97) 

and 

L = L' - Md^ - k^Ddl (6.98) 

L' = dt + L' + V^D. (6.99) 

Then (|3^ )-(|3l7|) become 

L'Uo,e = k:'Dd^eUo,e, (6.100) 

L'Uo,ee = F" : Ul^ + eDd',Uo,ee, (6.101) 

L'Uo,k = MUo,e + k^Dd''sUo,k, (6.102) 

L'Uo^kk = F":Ul, + 2MUo,ek + 2DUo,ee 

+ uUo,e + k^Dd^sUo,kk, (6.103) 

L'Uo^ek = MUo,ee + F":Uo,eUo,k + k'Dd',Uo,ek. (6.104) 



Suppose that S,U, and U are expanded in the same forms as (|5.64| ), (|5.65| ) 



and (|5.66| ), then up to O(e^) perturbed secular fields Uj (j = 1, 2, 3) obey 



L'Uj -- 


= Md^U,^^, 


(6.105) 


U, - 


= P,{x, r^, t)Uo,g + P,-_i,,t/o,fc + Uj, 


(6.106) 
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where Pq = Ui = 0. 

For j = 1, ( lETOSD and (|6106| ) give 



which imphes 
and 



Pi,tUo,e - V^PiDUo,e = 0, 

P,,t = 0, V^Pi = 0, 
Pi = Pi,.x + VxPi • r^. 



To j = 2, (|6105|) and (|O06| ) give 

L'f/2 = -P2,tt/o,e + \/^P2DUo,e. 
The similar discussion as in the previous sections yields 

U^ = S2W{e), 52 = V^P2 = constant, P2,i = 0, 
and 

P2 = P2,..xy2 + x{r^ : Vx)P2,.- + (r^r^ : Va_V±)P2/2, 

where P2,xx, ^±P2,x and (V_l)^P2 are arbitrary constants. 
Similarly, to j = 3, we have 

U3 = SsW{e), Ss = V^Ps = S3,xX + V±S3-r^, Ps.t = 0, (6.113) 



(6.107) 

(6.108) 
(6.109) 

(6.110) 
(6.111) 
(6.112) 



and 



P3,xxxxy3\ + x''{r^ : Vx)P3,../2 + xir^r^ : ViVx)P3,./2 



+ {r^r^r^ : V±V±V±)P3/3\, 



(6.114) 



which is a polynomial of degree three and all the coefficients are arbitrary 

constants 

To 0(e®), a nonlinear term appears as 

L't/4 = (1/2)P" : U^ + MdxU^, (6.115) 

t/4 = P4Uo,e + P3,xUo,k + il/2)P^UoM + U4i0,x, r^,t). (6.116) 

Substituting ( |6.116| ) into ( |6.115| ), we get 



+V^{P^/2)DUo,ee-Pi,tUo 



^o,e- 



Let set 



U^i9) = SiW{9) + Uiie), 



(6.117) 
(6.118) 
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where 

^4 = ^4,....xV2 + x(rx : Vx)^4,. + (r^r-x : Vi_V±)S^/2, (6.119) 

which is a polynomial of degee two, then we have 

L'U4 = -{S^-V^P4)DUo,e + V^Ps,,DUo,k + P2,xxMUo,k 

-V^S4DW - S3^^MW + \/'^{P^/2)DUo,e - P4,tUo,e- (6.120) 

This equation is valid only when all the coefficients of functions of 6 in ( |6.120|) 
are constant and the compatibility condition gives 

(V'P2/e')(f> • DUo,e) + P2,xx{U ■ MUo,k) + V^P^AU ■ MW) 
+'^^P3,x{U ■ DUo,k) + \VPif{U ■ DUo,ee) + V^^P^iU ■ DW) 
-P4AU-Uo,e) = 0. (6.121) 

Now a renormalization transformation to O(e^) takes the same form as ( |5.86| ). 
From (|5:86|) , i ^AWj ) , i ^JT% ) , ( |6Tl4D and ( |6A|) we obtain 



e^Pi,x, Vx0 = e^VxPi, 4>xx 


— C P2,xx, 


Vi0 = eViP2, 




= e^Ps,..., Vi4 = e^ViPa,., 


[^ tJL/ JU JU JU 


e Pa,xxxx-, V_|_0 = 


e'ViP4, 


e'P4,t. 






(6.122) 



Substituting ( |6.122[ ) into ( |6.121| ), we arrive at a generalized Benney equation 



in the three dimensional space: 

0t = I)|i0xx + /^±Vi0 + ^(0,,, + Vl4)x) + ^V"0 + iV±|V0p. (6.123) 

As in the case of the n-K-S equation, secular terms in U are also automatically 
removed when 6 is renormalized. 



7 Complex Ginzburg-Landau Equation 



As an application of the previous results, let us calculate explicitly various 
coefficients of the phase equations for the complex Ginzburg-Landau (cGL) 
equation. 
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where \1/ is a complex variable, 7 is a real constant, a and /3 are complex 
constants. The bar denotes complex conjugation. Setting 

the cGL equation is transformed into the standard form of a reaction-diffusion 
system ( |3.1| ). As a periodically oscillating solution of ( |7.124|) , we take 

f/o=(|)=a(A:)(^';l), (7.126) 

where a real amptitude a{k) and a frequency ooik) satisfy the following disper- 
sion relation 

- 2W = 7 - /3a^ - aA;2, (7.127) 

which Loik) and aik) are functions of k^ as speculated in section 6. Since 
the ajoint operator V' of L defined in (p.3|) is 

Lt = ujde - F'^iUo) - D^k^dl (7.128) 

where 

we explicitly find the null vector of L^ as 

U = t/2(3'a(^J^''^l,y (7.130) 

with the normalization {U ■ Uq^) = 1 and (3 = (3' + ij3" . Various differentials 
of Uq are also explicitly given by 

u I '^'' \ u ■ ( ^'' \ u I ^'' \ 

Uo,e - a I _-^-i0 I , tyo,fc - a I ^_jg 1 , Uo,e0 - -a I ^-ie I , 

Uo,ek = a i _^^^_ig j , f/o,fcfc = a i ^_ig j . (7.131) 

Using ( [7. 1301 ) and ( |7.131D , the coefficients of diffusion and nonlinearity are 
calculated as 

D^_ = {U-DUo,e) = Re{ap)/(3', (7.132) 

D|l = {U ■ MUo,k) = {-d/(3'a){dj(3" - 2kRe{a^)} 

= -2d2|/?|7/?', (7.133) 

N^ = (U . DUo,ee) = -Im{a^)/(3' = -u/2 = A^,,, (7.134) 
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where iVy happens to be identical with A^_l in this case. The coefficient of dis- 
persion in the K-dV-Burgers equation is calculated as follows. The assumption 
(g) reads 

D\\ = {Re{a(3) - 2a^\(3\^] / (3' ~ 0(e), (7.135) 

and 

i?e(a/3) ~2a2|/3|l (7.136) 

Then V in (|4.34|) is obtained as 

V{d) = {Im{a(3)/2ad\f3\'^)Uo,k, (7.137) 

and we have the coefficient of dispersion 

A' = {U-DUo,k) = aIm{a^)/a(3\ (7.138) 

B = {U- MV{e)) = {Im{a(3)/2ad\f3\'^){U ■ MUo,k) 

= -dIm{aP)/aP', (7.139) 

A = B + A' = 2dlm{a^)/a(3'. (7.140) 

From the Kuramoto-Sivashinsky scale ( p.61[ ), we have 

D^= {U ■ DUo,e) = Re{ap)/(3' ~ Re{ap) ~ 0{e^). (7.141) 

Then W in ( |5.63|) i s found to have the same expression as V since the forcing 
terms in Eqs.( |4.34D a nd ( p. 63 ) have the same non-seculer component under 
the conditions ( [7. 1351) and 



The coefficients of the K-S equation are obtained as 

E = {U- DW{e)) = {Im{a(3)/2ad\(3\'^){U ■ DUo,k) 

= -{Im{a(3))y2a^(3'\(3\\ (7.142) 

H = {U ■MW{9))=B, (7.143) 

G = A' + H = A' + B = A. (7.144) 

All the coefficients of the Benney equation are given in terms of those of the 
K-B equation and the K-S equation. The isotropic part of coefficients of the 
n-Burgers equation and the K-S equation agrees with the result in M. 



8 Concluding Remarks 



Let us compare the present derivation of phase equations with a possible 
derivation by means of the reductive perturbation (RP) method 0| or the 
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multi space-time scale method, although the latter derivation has not been 
accomplished yet. The initial setting of perturbation ( |3.8| ) (and ( p.9| )) and 
auxiliary conditions ( [4.32| ) and ( |5.61| ) are same for both derivations. In addi- 



tion to the initial setting and the auxiliary condition, the RG method assumes 
a naive expansion of secular solutions such as (|3.1CI|) - (|3.12|) and ( |4.35|) -( p:.37|) 



etc. . Then straightforward calculations of secular terms with the aid of the 
consistent renormalization condition and the RG procedure of the Lie approach 
lead to the final results. The type of derived phase equations depends only 
on a specific form of expansion of secular terms and the auxiliary condition. 
On the other hand, the RP method would requires specific scalings for not 
only perturbed fields but also independent variables, which are available after 
a derived equation is set up. It should be mentioned that the present RG 
method relies on an explicit secular solution although the RP method would 
not. However, the step to obtain an explicit secular solution may be largely 
skipped as far as the final results are concerned or possibly by the proto-RG 
approach developed recently M. 
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